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ABSTRACT 

Astrophysical black hole candidates are thought to be the Kerr black holes 
predicted by General Relativity, but the actual nature of these objects has still 
to be proven. The Kerr black hole hypothesis can be tested by observing strong 
gravity features and check if they are in agreement with the predictions of Gen- 
eral Relativity. In particular, the study of the properties of the electromagnetic 
radiation emitted by the gas of the accretion disk can provide information on the 
geometry of the space-time around these objects and constrain possible devia- 
tions from the Kerr background. 

Subject headings: black hole physics — general relativity — accretion, accretion 
disks 



Introduction 



The predictions of General Relativity have been confirmed by experiments in the Earth's 



gravitational field (Williams et al. 2004 Everitt et al. 2011), by spacecraft missions in the 



Solar System (Bertotti et al. 2003), and by accurate radio observations of binary pulsars 



(Weisberg & Taylor 2005 Kramer et al. 2006) (for a general review, see e.g. Will (2006)). 



In all these environments, the gravitational field is weak, in the sense that one can write 
Qtt = 1 + (p with 101 ^ 1. The validity of the theory in the regime of strong gravity, when 
the approximation |0| <^ 1 breaks down, is instead still unexplored. The ideal laboratory 
to test strong gravitational fields is the space-time around astrophysical black hole (BH) 
candidates. 
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In 4-dimensional General Relativity, uncharged BHs are described by the Kerr solution 
and are completely specified by two parameters: the mass, M, and the spin angular momen- 
tum, J (Carter 1971 Robinson 1975)^ A fundamental limit for a Kerr BH is the bound 
I a* I < I7 where a* = J/M^ is the spin parameter. This is the condition for the existence 
of the event horizon: for \a^\ > 1, the event horizon disappears, and the central singularity 
becomes naked, violating the weak cosmic censorship conjecture ( Penrose|1969 ). Despite the 



apparent possibility of forming naked singularities from regular initial data (see e.g. Joshi 
& Malafarina ( 2011a||b ) and references therein), the existence of a Kerr naked singularity or 
of a super-spinning Kerr compact object (in which the singularity may be replaced by an 
extremely compact object made of exotic matter) can be excluded for at least two reasons: 
it is apparently impossible to make a star collapse with |a^,| > 1 (Giacomazzo et al. 2011) or 



overspin an already existing Kerr BH up to |a=i,| > 1 (Barausse et al. 2010, 2011), and, even 



if created, this super-spinning object would be unstable and should quickly decay (Dotti et 



al. 


2008 


Pani et al. 


2010 



Astrophysical BHs presumably form from the gravitational collapse of matter, while the 
BH uniqueness theorems require a stationary space-time. However, initial deviations from 
the Kerr solution should be quickly radiated away through the emission of gravitational waves 
(Price 1972a||b ) (but see Lyutikov &; McKinney (2011)), and the geometry of the space-time 
around the compact object should be well described by the Kerr metric]^ At the moment. 



there are at least two classes of astrophysical BH candidates (for a review, see e.g. Narayan 
(2005)): stellar-mass compact objects in X-ray binary systems (M ^5 — 20 Mq 

10^ - 10^ Mr, 



and super- 
The measurement 



massive bodies at the center of every normal galaxy (M ~ ±u — ±u iviq 
of the mass of the two classes of BH candidates is robust, because obtained by dynamical 
methods and without any assumption about the nature of these objects: basically, we can 
study the orbital motion of gas or of individual stars orbiting around the compact object 
and we can infer the mass of the latter by using Newtonian mechanics. Stellar-mass BH 
candidates are too heavy to be neutron or quark stars ( Rhoades fc RufEiiI||1974 Kalogera & 
Baym|1996 ), while at least some of the super-massive objects in galactic nuclei are too heavy. 



compact, and old to be clusters of non-luminous bodies (Maoz 1998). The non-observation 



of electromagnetic radiation emitted by the possible surface of these objects may also be 



^The Kerr metric is a vacuum exact solution even in other theories of gravity, but its uniqueness is not 
guaranteed in general ( Psaltis et al.|200"8 ). The scenario of Kerr BHs in an alternative theory of gravity may 
be tested by looking for observational features that do not depend only on the metric of the background, 
such as the quasi- normal modes of the compact object ( Barausse fc Sotiriou||2008 ). 

^Deviations from the Kerr metric due to the presence of an accretion disk should also be extremely small 
in the case of a compact object in an X-ray binary system. However, the presence of surrounding material 
might be important for some super-massive compact objects in galactic nuclei. 
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interpreted as an indications that they have an event horizon and are therefore BHs (Narayan 



fc Heyl|2002]|McChntock et aL|2004[|Narayan fc McChntock|2008j [Broderick et al.|2009| ) (see 



however Abramowicz et aL (2002) and Bambi (2012e)). As there is no alternative explanation 
in the framework of conventional physics, astrophysical BH candidates are thought to be the 
Kerr BHs predicted by General Relativity, even if we have no evidence that the geometry of 
the space-time around them is really described by the Kerr solution. 

The first proposal for testing the Kerr BH hypothesis was put forward by Ryan, who sug- 
gested to observe the gravitational waves emitted by an extreme-mass ratio inspiral (EMRI), 
i.e. a system consisting of a stellar-mass compact object orbiting around a super-massive 



BH candidate (Ryan 1995): as future space-based gravitational wave detectors will be able 
to observe ~ 10^ — 10^ gravitational wave cycles emitted by an EMRI while the stellar-mass 
body is inspiraling into the gravitational field of the super-massive object, even a small de- 
viation from the Kerr geometry will build up an observable dephasing in the gravitational 
waveforms, thus allowing one to map the space-time of the super-massive BH candidate 



with very high accuracy (Glampedakis & Babak 2006 Barack &; Cutler 2007). Besides tests 



based on gravitational waves, the space-time geometry around BH candidates can be probed 
by very accurate observations of the orbital motion of stars. In the case of a stellar-mass 
BH candidate in an X-ray binary system, that is possible if the companion star is a radio 
pulsar ( Wex fc Kopeikin||1999 ). In the case of the super-massive BH candidate at the center 
of our Galaxy, that might be achieved by monitoring stars orbiting at milliparsec distances 



from the compact object (Will 2008). 



Recently, there has been an increasing new interest in the subject. While previous stud- 
ies focused on future facilities (like space-based gravitational wave detectors) or possible fu- 
ture discoveries (like a BH binary with a pulsar companion), more recent efforts deal with the 
possibility of testing the Kerr-nature of astrophysical BH candidates by studying the proper- 
ties of the electromagnetic radiation emitted by the gas of the accretion disk, for which there 
are already available data. In the last 5-8 years, there have been significant progresses in the 
understanding of the accretion process onto BH candidates and astronomers have developed 
several approaches to estimate the spin parameter of these objects under the assumption of 
Kerr background. One can typically extend these techniques to constrain/discover possible 



deviations from the Kerr metric (Bambi 2011d). 



The content of the paper is as follows. In Section |2| I review the theoretical framework 
to test the Kerr-nature of astrophysical BH candidates: the Novikov-Thorne model, which 
can be used to describe thin accretion disks in generic stationary and axisymmetric space- 
times, and the choice of a non-Kerr background metric to parametrize possible deviations 
from the Kerr solution. In Section |3| I show some approaches to probe the geometry of 
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the space-time around astrophysical BH candidates from the analysis of the electromagnetic 
radiation emitted by the gas of accretion: the estimate of the radiative efficiency of AGNs, 
which can provide a rough idea of the size of deformations that may already be considered 
unlikely, the continuum-fitting method and the broad Ka iron line analysis, which are the 
two most popular techniques currently used by astronomers to estimate the spin parameter 
of these objects and can be naturally extended to constrain possible deviations from the Kerr 
background, the high-frequency quasi-periodic oscillations (QPOs) and the estimate of the jet 
power, which may potentially be powerful probes, but in both cases the associated physical 
mechanism is not yet well understood from the theoretical point of view. Throughout the 
paper, I use units in which Gn = c = 1, unless stated otherwise. 



2. Accretion process onto a compact object 



2.1. Novikov-Thorne model 



Geometrically thin and optically thick accretion disks are described by the Novikov- 
Thorne (NT) model ( Novikov fc Thoriie||1973 Page fc Thorne||1974 ), which is the relativistic 



generalization of the Shakura-Sunyaev model (Shakura & Sunyaev 1973). The disk is thin 
in the sense that the disk opening angle is h = H/r <C 1, where H is the thickness of 
the disk at the radius r. In the Kerr background, there are four parameters (BH mass 
M, BH spin parameter a^,, mass accretion rate M, and viscosity parameter a), but the 
model can be easily extended to any (quasi-) stationary, axisymmetric, and asymptotically 
fiat space-time. Accretion is possible because viscous magnetic/turbulent stresses and ra- 
diation transport energy and angular momentum outwards. The model assumes that the 
disk is on the equatorial plane and that the disk's gas moves on nearly geodesic circular 
orbits. For long-term accretions, the disk is adjusted on the equatorial plane as a result 



of the Bardeen-Petterson effect (Bardeen & Petterson 1975) (such a mechanism works even 



in non-Kerr backgrounds (Bambi 2011a)). In the case of BH candidates in X-ray binaries. 



the disk is likely on the equatorial plane if the two bodies formed from the collapse of the 



same cloud (Fragos et al. 2010). The assumption of nearly geodesic circular orbits requires 



that the radial pressure is negligible compared to the gravitational force of the BH, which 
is indeed the case for thin disks. Heat advection is ignored (it scales as ~ h'^) and energy is 
radiated from the disk surface. Magnetic fields are also ignored. 

The key-ingredient of the NT model is that the inner edge of the disk is at the innermost 
stable circular orbit (ISCO), where viscous stresses are assumed to vanish. When the gas's 
particles reach the ISCO, they quickly plunge into the BH, without emitting additional 
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radiation. At first approximation, the total efficiency of tlie accretion process is 

r^tot = 1 - ^isco , (1) 

wfiere -Eisco is tlie specific energy of ttie gas at tlie ISCO radius and depends uniquely on tlie 
background geometry (see Appendix). In general, tlie total power of the accretion process 
is converted into radiation and kinetic energy of jet/wind outflows, so we can write ?7tot = 
^rad + ^kin, but in the NT model r^km = 0. ?7rad is the radiative efficiency and can be inferred 
from the bolometric luminosity Lboi if the mass accretion rate is known: Lboi = VradM. 

How good is the NT model to describe the accretion disk around astrophysical BH 
candidates? For non-magnetized and weakly-magnetized disks, there is a common consensus 
that the NT model describes correctly thin disks, /i <^ 1, when the viscosity parameter is 



small, a -C 1 (Afshordi & Paczyhski 2003). A common criterion to select sources with thin 



disks is that the bolometric to Eddington luminosity ratio, LboiZ-^^Edd, does not exceed 0.3 
dMcClintock et al.|[2006| . In the case of magnetized disks, the issue is open and controversial. 
The GRMHD simulations in Penna et al. (2010) (see also Penna et al. ( 2012[ )) show that the 
stress at the inner edge of the disk scales as h and the authors conclude that the NT model 
with vanishing stress boundary condition is recovered for — )■ 0. The GRMHD simulations 
in Noble & Krolik (2009) and Noble et al. (2010) show instead large stress at the inner edge 



of the disk even when h ^ 0; according to these authors, the NT model cannot describe 
magnetized disks, even when the disk is very thin. It is not clear if the disagreement between 
the two groups can be attributed to different configurations of the magnetic fields, different 
resolution of the simulations, or something else. 



2.2. Parametrizing deviations from the Kerr geometry 

In the weak-field limit. General Relativity has been successfully tested in many sit- 
uations. Solar System's experiments are conveniently discussed within the parametrized 
post-Newtonian (PPN) formalism. One can see that the weak-field limit of a general static 
and spherically symmetric metric can be put in the form 

ds^ = -(l- — + 2(3^ + ..)dt^+(l + 2^— + ..)idx^ + dy^ + dz^), (2) 



where /3 and 7 are the PPN parameters. In General Relativity, /3 = 7 = 1, but in other 
metric theories of gravity that may not be true. Observations can measure (3 and 7 and thus 
constrain possible deviations from the predictions of General Relativity. 



The natural extension of the PPN formalism to test the space-time geometry around 
astrophysical BH candidates, where gravity is strong, is to write a general stationary, ax- 
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isymmetric, and asymptotically flat metric in terms of the mass moments, Aii, and of the 



current moments, Si, of the compact object (Ryan 1995). In the case of a Kerr BH, all the 



moments are locked to the mass and angular momentum by the following relation (Hansen 



1974): 



Ml + iSi = M 



iJ_ 
M 



(3) 



By measuring at least three non-trivial multiple moments of the space-time around an as- 
trophysical BH candidate, one can test the Kerr-nature of the body. The advantage of a 
post-Newtonian approach is that we do not need specific assumptions about the nature of 
the BH candidate. The disadvantage is that we can only probe the metric relatively far 
from the compact object, where the contribution of higher-order moments can be neglected, 
and therefore we need very accurate measurements, which may not be possible in a realistic 
astrophysical situation. 

In order to test the space-time geometry close to a BH candidate, we can consider a 
metric describing the gravitational field (in General Relativity or in an alternative theory 
of gravity) of a compact object characterized by its mass, its spin angular momentum, and 
a set of deformation parameters measuring possible deviations from the Kerr solution. In 
this way, one can compare observations and theoretical predictions in this more general 
background and check if the observations require that all the deformation parameters must 
vanish, as necessary if the space-time geometry around the BH candidate is really described 
by the Kerr metric. In the simplest case, we can consider just one deformation parameter, 
which should be able to describe the gravitational field around a compact object either 
more oblate or more prolate (according to the sign of the deformation parameter) than a 
Kerr BH. If our metric does it, it is not very important if it is an exact solution of the 
vacuum Einstein's equations (in this case, we are actually considering the possibility the BH 
candidate is a compact object made of exotic matter) or not, in the sense that we eventually 
get the same qualitative constraints ( Bambi|[2011e 2012a Li &: Bambi||2012 ). In this review 
paper, I will focus the discussion on the metric proposed by Johannsen and Psaltis (JP) 
to describe the gravitational field around non-Kerr BHs in putative alternative theories of 
gravity ( Johannsen &: Psaltis||2011b ). In Boyer-Lindquist coordinates, the JP metric is given 
by the line element 



f 2Mr\ , , , . 
(^1-— j (l + /i)rft2 + 



S(l + /i) 



+ 



A + a^h sin^ Q 
sin^^ 



dr^ + S dQ'^ 



AaMr sin^ 6 



l + h)dtd(f) + 



2 I 2 

r + a 



2a^Mr sin^ 9 



+ 



(S + 2Mr) sin^ d 



h 



(4) 
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where a = a^,M, S = + cos^ 6, A 



h 



E 

fc=0 



— 2Mr + a^, and 
Mr 



e2fc 



M2 

IT 



(5) 



This metric has an infinite number of deformation parameters and the Kerr solution 
is recovered when all the deformation parameters are set to zero. However, in order to 
reproduce the correct Newtonian limit, we have to impose eo = ei = 0, while €2 is strongly 
constrained by Solar System experiments ( Johannsen & Psaltis 2011b ). In what follows, I will 
examine only the simplest case in which €3 7^ 0, while all the other deformation parameters 
are set to zero. 



3. Testing the Kerr-nature of BH candidates 
3.1. Radiative efficiency of AGNs 



The value of the spin parameter of a compact object is determined by the competition of 
three physical processes: the event creating the object, mergers, and gas accretion. Accretion 
from a disk can potentially be a very efficient way to spin a compact object up. If the inner 
edge of the disk is at the ISCO radius, the gas's particles plunge into the compact object 
with specific energy -Eisco and specific angular momentum Lisco- The mass M and the spin 
angular momentum J of the compact object change respectively by SM = EjscoSfn and 
6 J = Lisco^m, where Sm is the gas rest-mass. The evolution of the spin parameter turns 



out to be governed by the following equation (Bardeen 1970) 

da^ 1 L 



ISCO 



dlnM M 



2a* 



ISCO 



(6) 



neglecting the small effect of the radiation emitted by the disk and captured by the object. 
Prolonged disk accretion is a very efficient mechanism to spin the compact object up till 
an equilibrium spin parameter a^'^, which is reached when the right-hand side of Eq. (|6]) 
vanishes. For instance, an initially non-rotating Kerr BH reaches the equilibrium a^^ = 1 
after having increased its mass by a factor ^ 2.4 teardeen||l970|). 



One can argue that the most favorable scenario to produce fast-rotating super-massive 
objects at the center of galaxies is likely via prolonged disk accretion and that the maximum 



value for the spin parameter of these objects today cannot exceed (Bambi 201 le). The 
numerical value of a^^ depends on the metric of the space-time and it may exceed 1, as the 



bound I a* I < 1 holds for Kerr BHs ( Bambi|2011a|[c Bambi &: Modesto|2011 ). The black solid 
curve in Figure [T] shows the equilibrium spin parameter for the JP metric with deformation 
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parameter e^. Objects on the left of the black solid curve have a* < al^ and the accretion 
process spins them up; objects on the right have a^, > af^ and the accretion process spins 
them down. 

In general, it is not easy to get an estimate of the radiative efficiency r^rad, as the 
measurement of the mass accretion rate M is typically quite problematic. The mean radiative 



efficiency of AGNs can be inferred from the Soltan's argument (Soltan 1982), which relates 



the mass density of the super-massive BH candidates in the contemporary Universe with the 
energy density of the radiation produced in the whole history of the Universe by the accretion 
process onto these objects. There are several sources of uncertainty in the final result, but 



a mean radiative efficiency ?7i.ad > 0.15 seems to be a conservative lower limit (Elvis et al. 
2002). Wang et al. (2006) find a mean radiative efficiency r/rad ~ 0.30 — 0.35 without some 



important assumptions necessary in the original version of the Soltan's argument. Davis fc 



Laor (2011) show how to estimate ?7rad for individual AGNs and find that the most massive 
objects have typically higher ?7rad, up to ~ 0.3 — 0.4. 

It seems plausible to conclude that at least some of the super-massive BH candidates 
must have r^rad > 0.15. While the accretion disk of these objects may not be necessary 
thin, accretion from a thin disk is among the most efficient ways to convert the rest-mass 
of the gas into radiation. As r/tot > ''7rad, -Eisco < 0.85 and we can constrain deviations 



from the Kerr geometry (Bambi 2011b|e). A simplified explanation of how this is possible is 



the following. Ejsco and a^^ increase/decrease if the gravitational force around the object 
increases/decreases, because the latter moves the ISCO to larger/ smaller radii. Under the 
conservative assumption that ?7tot = Vrad, we can compute r^tot of the background metric 
and, combining this bound with the one for the maximum value of the spin parameter 
as a function of the deformation parameter, we get a well defined allowed region in the 
spin parameter-deformation parameter plane. In the specific case of the JP metric with 
deformation parameter es, we find the plot in Figure [TJ While the specific value of the 
deformation parameter depends on the choice of the metric of the background and a physical 
interpretation is not straightforward, the maximum value for the spin parameter seems to 
depend only very weakly on the exact space-time, which makes this argument quite appealing 



(Bambi 2011e). The maximum value today of the spin parameter for the super-massive BH 



candidates in galactic nuclei would be 1.19 if we require T^tot > 0.15, 1.10 for ijtot > 0.20, and 
1.04 for r/tot > 0.25. 



More recently, Li & Bambi (2012) have studied the accretion process from thick disks. 



While in this case the details strongly depend on the unknown properties of the accretion flow, 
at least in principle a thick disk can spin up a BH candidate more efficiently. The evolution 
of the spin parameter is still given by Eq. ([6]), but the quantity i^isco/-Eisco is replaced 
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by Li^J Ei^^, where i?in and Ljn are, respectively, the specific energy and the specific angular 
momentum of the gas particles at the inner edge of the disk. For a thick disk, the inner edge 
of the disk can be inside the ISCO. So, if the super-massive BH candidates have experienced a 
period of super-Eddington accretion in the recent past, their spin parameter might exceed the 
equilibrium value a^'^ computed from a thin disk. Under the most conservative assumptions, 
the previous bounds on the maximum value of become shghtly weaker: 1.30 for r^tot > 0.15, 
1.18 for ritot > 0.20, and 1.10 for r/tot > 0.25 (iLi & Bambi||2012[). 



3.2. Continuum-fitting method 



In the NT model, the thermal spectrum of a geometrically thin and optically thick accre- 
tion disk around a Kerr BH observed far from the compact object depends on 5 parameters: 
BH mass M, BH distance d, inclination angle of the disk with respect to the line of sight 
of the observer i, BH mass accretion rate M, and BH spin a^,. If we can get independent 
measurements of M, d, and i, we can fit the data and estimate M and a*. The technique is 



called continuum-fitting method (Zhang et al. 1997 Li et al. 2005; McClintock et al. 2011) 



and it can be applied only to stellar-mass BH candidates: the disk's temperature goes like 
M~°'^^, and the peak of the spectrum is at about 1 keV for M ~ 10 Mq and in the UV 
range for super-massive objects in galactic nuclei. In the latter case, a good measurement is 
impossible because of dust absorption. From the conservation laws for the rest-mass, angular 
momentum and energy, one can deduce three basic equations for the time-averaged radial 
structure of the disk. In particular, the time-averaged energy flux emitted from the surface 



of the disk is (Page & Thorne 1974) 



M 



47rM2 

where F{r) is the dimensionless function 



F{r) 



(7) 



F{r) 



{E - my 



{E ~nL){dpL)dp. 



(8) 



E, L, and Q are, respectively, the conserved specific energy, the conserved axial-component 
of the specific angular momentum, and the angular velocity for equatorial circular geodesies 
(see Appendix); G = —a'^grrdijxj) is the determinant of the near equatorial plane metric, where 



— 9t<f>/ ~ 9tt is lapse function; r^^ is the inner radius of the accretion disk and in 
the NT model it is the ISCO radius. 
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Since the disk is in thermal equihbrium, the emission is blackbody-hke and we can 
define an effective temperature Teff(r) from the relation J^{r) = aTjj, where cr is the Stefan- 
Boltzmann constant. Actually, the disk's temperature near the inner edge of the disk can 
be high, up to ~ 10^ K for stellar-mass BH candidates, and non-thermal effects are non- 
negligible. That is usually taken into account by introducing the color factor (or hardening 
factor) /col. The color temperature is Tcoi(r) = fcoiTes and the local specific intensity of the 
radiation emitted by the disk is 



/e(^e) 



2E^ 1 



T 



^ ^colexp 



(9) 



where Ee is the photon energy, h is the Planck's constant, c is the speed of light, /cb is 
the Boltzmann constant, and T is a function of the angle between the 3- momentum of the 
photon emitted by the disk and the normal of the disk surface, say ^. The two most common 
options are T = 1 (isotropic emission) and = -k + j cos^ (limb-darkened emission). 



The calculation of the thermal spectrum of a thin accretion disk has been extensively 



discussed in the literature; see e.g. Li et al. (2005) and references therein for the Kerr case 



and Bambi & Barausse (2011a), Chen & Jing (2012a) and Bambi (20121) for a background 



with generic deviations from the Kerr solution. The spectrum can be conveniently written 
in terms of the photon fiux number density as measured by a distant observer, Ne^^^^- In 
the case of non-Kerr background, it is convenient to use a ray-tracing approach. The initial 
conditions (to, ''"o, ^o, 4>o) for the photon with Cartesian coordinates (X, Y) on the image plane 



of the distant observer are given by (Johannsen & Psaltis 2010) 



to 
ro 

Oo 



0, 

Y sin i + d cos i 

VX^ +Y^ + d^ 
X 

d sin i — Y cos i 



arccos 



arctan 



(10) 
(11) 
(12) 

(13) 



As the initial 3-momentum ko must be perpendicular to the plane of the image of the 
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observer, the initial conditions for the 4-momentum of the photon are 

K = ^ Ikol, (14) 

nna 9 — fl Ysini+dcosi 

< = / '^:^^^^ =^^|ko|, (15) 



a/X^ + (d sin i — y cos i)^ 
Xsini 
+ (dsini — y cosi)^ 



= , ..l""',.,.... !".!. (16) 



where /cg is inferred from the condition g^j,uk^k^ = with the metric tensor of a flat space-time 
(as the observer is located far from the compact object). The photon trajectory is numerically 
integrated backwards in time to the point of the photon emission on the accretion disk: in 
this way, we get the radial coordinate at which the photon was emitted and the angle ^ 
between the 3-momentum of the photon and the normal of the disk surface (necessary to 
compute T). 

The observer's sky is divided into a number of small elements and the ray-tracing pro- 
cedure provides the observed flux density from each element; summing up all the elements, 
we get the total observed flux density of the disk. In the case of Kerr background, one 
can actually exploit the special properties of the Kerr solution and solve a simplified set of 
differential equations. That is not possible in a generic non-Kerr background and one has to 
solve the second-order photon geodesic equations of the space-time. The photon flux number 
density is given by 



-^^^obs — ~^ / -^obs(-E'obs)c^^^obs = / 9 h{Ee)dQ,obs 

-C/obs J -c/obs 



A 



^ VkeV J J M 



M2 

exp 



TdXdY , , 

(18) 



( i^obs '\ 

V keV / 



- 1 



where /obs and £^obs are, respectively, the specific intensity of the radiation and the photon 
energy measured by the distant observer. dQohs = dXdY/ d"^ is the element of the solid angle 
subtended by the image of the disk on the observer's sky. g is the redshift factor 

^ _ -E'obs _ ^g'^obs ^j^g-j 
Ee kpUe 

where A;" is the 4-momentum of the photon, u'^^^ = (—1,0,0,0) is the 4- velocity of the distant 
observer, and — {ul, 0, 0, flul) is the 4- velocity of the emitter. Ie{Ee)/E^ = Iobs{Eobs) / E^^^^ 
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follows from the Liouville's theorem. Ai and A2 are given by (for the sake of clarity, here I 
show explicitly Gn and c) 



2{keVf fGnM 



0.07205 ( mV /kpc^ ^ 



7 ke V cm ^ s ^ 



\M^) \ d 
keV \ /GnM\^/^ fA-na^ 



Ao 



/col \ M J \Mq 

Using the normalization condition g^^u'^u'^ = — 1, one finds 



^b/coi/ \ J \ M 
0.1331 fW^gs-^V'^ { M^^'^ 



(20) 



< = , ^ , (21) 

\/-gtt - '^9tS - dH^"^ 

and therefore 

^~ i + xn ' ^^^^ 

where X = k^/kt is a constant of the motion along the photon path and can be evaluated 
from the photon initial conditions {X = k^/kt = tq sin^o^o/^o)- 

Figure [2] shows the thermal spectrum of a thin disk as a function of the parameters of 
the model. M, d, and i should be inferred from independent measurements, while a*, M, and 
es can be obtained by fitting the disk's spectrum. It is evident that variations of a* and €3 
produce quite similar effects on the shape of the spectrum, while M changes the spectrum in 
a very different way. Figure |3] shows the effects of variations of the color factor f^oi (typically 
one should expect /coi = 1-6 — 1.7) and of the function T (whose effect is relevant only for 
large inclination angles). 

In order to figure out how the analysis of the thermal spectrum of a thin disk can 
constrain the geometry of the space-time around a BH candidate, we can compare the 
spectra of a Kerr BH and of a JP BH with deformation parameter e^. We define the reduced 
as follows 

2 



2 1 - iV/P(a.,e3,M)-iVf-'^(S„M) 
xUa*, ^3, M) = ^ = - 5^ i L , (23) 



where the summation is performed over n sampling energies Ei, N^^ and N^'^'^^ are the photon 
flux number densities respectively for the JP and Kerr metric, and (jj is the uncertainty. As 
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variations of M produce quite different effects on the shape of the spectrum with respect to 
a^, and £3, its determination is not correlated to the other two parameters and here we directly 
assume M = M. Figure [4] shows xicd ^"^^ cases: a^, = 0.7 (left panel) and a* = 0.98 
(right panel), assuming that the uncertainty is 10% the photon flux number density; that 
is, (jj = 0.1 Nf^'^". The continuum-fitting method basically measures the radiative efficency 
^rad = 1 — Eisco, not the Spin parameter a^,. In the Kerr metric, ?7rad depends only on a^: and 
therefore one can infer the spin of the compact object. In a generic non-Kerr background, one 
can only measure a combination of a^, and of the deformation parameters. The technique 
is indeed unable to distinguish a Kerr BH with spin parameter a^, from a very deformed 



object with completely different spin parameter but very similar radiative efficiency (Bambi 



fc Barausse||2011a ). An analysis of real data from specific sources will be presented in Bambi 



et al. (2013). 



For the time being, the continuum-fitting method has been used to estimate the spin 
parameter of 9 stellar-mass BH candidates under the assumption that these objects are 
Kerr BHs. The list of these objects is reported in Table [T| There are 2 objects that look 
like very fast-rotating Kerr BHs (GRS 1915-1-105 and Cygnus X-1), 4 objects that seem to 
be Kerr BHs with a mid value of the spin parameter (LMC X-1, M33 X-7, 4U 1543-47, 
and GRO J1655-40), and 3 objects whose data are consistent with slow-rotating Kerr BHs 
(XTE J1550-564, LMC X-3, and A0620-00). The measure of the spin parameter a^, under 
the Kerr BH hypothesis can be quickly translated into a measurement of ?7rad and then 
into an allowed region on the spin parameter-deformation parameter plane, see Fig. [5j The 
comparison of the allowed region for GRO J1655-40 (black area, right panel in Figure [s]) 
with the Xred ^^^^ panel of Figure [4] shows that the continuum- fitting method measures 

indeed the radiative efficiency. For high values of the spin and of the deformation parameters, 
the rule of the radiative efficiency does not work so well (see the case of GRS 1915-1-105 or 
Cygnus X-1 and compare with the allowed region in the right panel of Figure |4]), but such a 
region of the plane is more likely unphysical, as it is probably impossible to spin these very 
prolate compact objects to so high values of a* (see e.g. the black line in Figure [T]). 



3.3. Broad Ka iron line 

The X-ray spectrum of both stellar-mass and super-massive BH candidates is often char- 
acterized by the presence of a power- law component. This feature is commonly interpreted 
as the inverse Compton scattering of thermal photons by electrons in a hot corona above 
the accretion disk. The geometry of the corona is not known and several models have been 
proposed. Such a "primary component" irradiates the accretion disk, producing a "reflec- 
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0.5 1 1.5 2 1 1.1 1.2 1.3 



a. a. 

Fig. 1. — Total efficiency r^tot = 1 — -E'isco in the JP background with spin parameter a^, and 
deformation parameter €3. r/tot = 0.15 (red solid curve), 0.20 (green dashed curve), and 0.25 
(blue dotted curve). The black solid curve is the equilibrium spin parameter al'^ obtained 
from Eq. (|6|. The right panel is an enlargement of the parameter region 1.0 < a^: < 1.3 and 
-2.0 < £3 < 0.0. 



BH Binary 



-.Kerr 



Vn 



Vn 



Reference 



GRS 1915+105 
Cygnus X-1 



a, > 0.98 
a, > 0.97 



0.234 
0.215 



0.423 
0.423 



McClintock et al. 



2006 





Gou et al. 


(2C 


)11 


) 




Gou et al. 


(2009 


) 


Li 


u et al. 


(20 


08 


20 


10) 



LMC X-1 
M33 X-7 
4U 1543-47 
GRO J 1655-40 



0.92 ±0.06 
0.84 ±0.05 
0.80 ±0.05 
0.70 ±0.05 



0.139 
0.120 
0.112 
0.097 



0.234 
0.151 
0.136 
0.112 



Shafee et al. 



Shafee et al. 



(2006) 



(2006) 



XTE J1550-564 
LMC X-3 
A0620-00 



0.34 ±0.24 

a* < 0.3 
0.12 ±0.19 



0.0606 
0.0365 
0.0550 



0.0892 
0.0694 
0.0699 



Steiner et al. 020111) 



Davis et al. 



Gou et al. 



(2006) 



2010 



Table 1: Continuum-fitting method results from the Harvard-Smithsonian CfA group. See 
references in the last column for more details. 
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log Eobs log Eobs 

Fig. 2. — Thermal spectrum of a thin disk as a function of the parameters of the modeL 
Left top paneh mass of the BH candidate in units of Mq. Right top paneh distance of the 
BH candidate in kpc. Left central panel: inclination angle of the disk with respect to the 
line of sight of the distant observer. Right central panel: mass accretion rate in units of 
]^gi8 g g-i bottom panel: spin parameter of the BH candidate. Right bottom panel: 

deformation parameter of the BH candidate. When not shown, the value of the parameters 
is: M = 10 Mq, = 10 kpc, i = 45°, M = 2- 10^^ g s-\ a, = 0.7, eg = 0.0. Flux density 
NE^^^^ in 7 keV~^ cm~^ s~^; photon energy Eohs in keV. 
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tion component" in the X-ray spectrum. The illumination of the cold disk by the primary 
component produces also spectral lines by fluorescence. The strongest line is the Ka iron 
line at 6.4 keV. This line is intrinsically narrow in frequency, while the one observed appears 
broadened and skewed. The interpretation is that the line is strongly altered by special and 
general relativistic effects, which produce a characteristic profile, first predicted in [Fabian 



et al. (1989) and then observed for the first time in the ASCA data of the Seyfert 1 galaxy 



MCG-6-30-15 (Tanaka et al. 1995). In MCG-6-30-15, this line is extraordinary stable, in 



spite of a substantial variability of the continuum, suggesting that the analysis of its shape 
can be used to probe the geometry of the space-time around the BH candidate. However, the 
relativistic origin of the observed broad Ka iron lines is not universally accepted, and some 



authors have proposed different explanations (Turner et al. 2007 Titarchuk et al. 2009) 



Within the interpretation of a relativistically broadened Ka iron line, the shape of 
the line is primarily determined by the background metric, the geometry of the emitting 
region, the disk emissivity, and the disk's inclination angle with respect to the line of sight 



of the distant observer (Fabian et al. 1989; Reynolds & Nowak 2003 Miller 2007). The 



effects of deviations from the Kerr metric on the shape of the Ka iron line are discussed 
Psaltis & Johannsenl ([20T2|, iJohannsen & Psaltisl (|20T2|, and iBambil (120131). In Kerr 



m 



and non-Kerr backgrounds, M sets the length of the system, so everything scales as M 
or as some power of M, without affecting the shape of the line. The relevant parameters 
of the background geometry are thus the spin a^: and possible non-vanishing deformation 
parameters. In those sources for which there is indication that the Ka iron line is mainly 
emitted close to the compact object, the emission region may be thought to range from 
the ISCO radius, nn = risco, to some outer radius rout- However, even more complicated 



geometries have been proposed (Martocchia et al. 2000). In principle, the disk emissivity 



may be theoretically calculated. In practice, that is not feasible at present. The simplest 
choice is an intensity profile Je cx r", with a < a free parameter to be determined during 
the fitting procedure. 

The line spectra can be expressed in terms of the photon flux number density as mea- 
sured by the distant observer: 

1 /■ . , „ s _ 1 



(-£'obs)o^^^obs 



(24) 



-^obs J -^obs 

Let US assume that the disk emission is monochromatic (the rest frame energy is i?Ka = 
6.4 keV) and isotropic with a power-law radial profile: 

Je(Se) oc 6iE, - EKa)r" . (25) 

The resulting broad Ka iron lines for different values of the model parameters (spin a^,, 
deformation €3, inclination angle i, power-law index a, and outer radius rout) are shown in 
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Figure |6} The photon flux has been normalized so that 

NE^,^JEohs = constant , (26) 



as only the shape matters. 

As we can see in Figure |6} the effect of variations of a* and of the deformation parameter 
€3 is similar. For small deformations, backgrounds with the same efficiency in the NT model, 
^rad = 1 — -E'iscO; have similar broad Ka iron lines and it is difficult to distinguish a Kerr 
BH from another compact object (see the right bottom panel). However, unlike in the case 
of the thermal spectrum of a thin disk, for larger deformations that is not true any more. 
The location of the ISCO radius (set by a* and €3 and determining r/rad) is encoded in the 
low-energy tail of the broad line. It consists of photons emitted from the receding part of the 
accretion disk close to the inner edge and thus strongly affected by both gravitational and 
Doppler redshift. The main peak of the line is produced by the effect of Doppler blueshift, 
relatively far from the compact object, where the effect of gravitational redshift is not so 
strong. For small values of the deformation parameters, the Keplerian velocity profile is still 
similar to the one of a Kerr BH and an object with similar radiative efficiency may look like a 
Kerr BH with different spin parameter. For larger deviations from the Kerr background, the 
role of the Doppler boosting (in both the main peak and the low-energy tail) becomes more 
and more important and breaks the degeneracy between spin and deformation parameter. 
To be more quantitative, we can proceed as in the previous section and compare the line 
produced in the space-time around a Kerr BH and a deformed object. The reduced is 
now: 



2 . . N 1 ^ [iV/P(a„ 63, i, a, Tout) - A^P"(a*,«, a, rout)] ^ 

Xrcd(a*> es, I, a, rout) = — = - > ^ 2 ' \^'> 

■ =1 



where the summation is performed over n sampling energies Ei, N^^ and Nf^"^" are the 
normalized photon fluxes in the energy bin [Ei, Ei + AE], respectively for the JP and Kerr 
metric, and cxj is the uncertainty of the energy bin. Motivated by the fact that variations 
of i, a, and rout produce different effects on the shape of the fine with respect to a^: and 
€3, for the sake of simplicity we can fix i = i, a = a, and rout = ''"out (actually, such an 
approximation may not be so good for the inclination angle, see Bambi| (2013)). Assuming 



a 10% uncertainty in the photon count, i.e. cr, = 0.1 iV^^*^", and data in the energy range 
1 — 7.5 keV with an energy resolution AE = 100 eV, we get the plots in Figure [7] for d = 0.7 
(left panel) and 0.98 (right panel). Unlike for the analysis of the disk's spectrum, here we can 
exclude that a source with a Ka iron line that looks like the one generated in a Kerr space- 
time with spin parameter a = 0.7 comes instead from the region around a very deformed 
compact object. 
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It is also interesting to see what happens if we combine the analysis of the thermal 
spectrum of the disk with the one of the broad Ka iron line. In general, it is not common 
to have good measurements of both the disk's spectrum and the Ka iron line for the same 
object, but it is possible and it may be more likely with future experiments. If we define 



Xred, 



tot 



2 _|_ 2 

Xred, Ka ' Xred, cfm ' 



(28) 



where Xred.Ka Xred,cfm respectively, the reduced defined in Eq. (27) and (23), we 
find the two plots shown in Figure |8| 



3.4. Quasi Periodic Oscillations 

The X-ray power density spectra of some low-mass X-ray binaries show some peaks; 



that is, there are QPOs in the X-ray flux (Remillard & McClintock 2006). These features 



can be observed in systems with both BH candidates and neutron stars and the frequencies 
of these oscillations are in the range 0.1 Hz-1 kHz. High-frequency QPOs in BH candi- 
dates (~40-450 Hz) are particularly interesting, as they depend only very weakly on the 
observed X-ray flux and for this reason it is thought they are determined by the metric of 
the space-time rather than by the properties of the accretion flow. If that is correct, they 
may be used to probe the geometry around stellar-mass BH candidates. For the time being, 
however, the exact physical mechanism responsible for the production of the high-frequency 
QPOs is not known and several different scenarios have been proposed, including hot-spot 



models (Stella & Vietri 1998 


199C 


)), diskoseismology models 


(Perez et al. 1997 Silbergleit 


et al. 


2001 Wagoner et al. 2001 ] 


"Cato 2001), resonance models (Abramowicz & Kluzniak 


2001 


Abramowicz et al.||2003 


Kluzniak & Abramowicz|2005 ; 


Torok et al.|2005), and p-mode 


oscillations of a small accretion torus (RezzoUa et al.||2003a|b 


Schnittman & Rezzolla||2006). 



In these models, the frequencies of the QPOs are directly related to the three characteristic 



BH binary 



M/M( 







z/u/Hz i^l/Hz 



References 



GRO J1655-40 
XTE J1550-564 
GRS 1915+105 



6.30 ±0.27 
9.1 ±0.6 
14.0 ±4.4 



450 ±3 
276 ±3 
168 ±3 



300 ±5 
184 ± 5 
113 ±5 



Strohmayer (2001) 



Remillard et al. 



(2002) 



Remillard & McClintock 



(2006) 



Table 2: Stellar- mass BH candidates in binary systems with a measurement of the mass and 
two observed high-frequency QPOs. 
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orbital frequencies of a test-particle: the Keplerian frequency z/k (which is the inverse of the 
orbital period), the radial epicyclic frequency z/,, (the frequency of radial oscillations around 
the mean orbit), and the vertical epicyclic frequency ug (the frequency of vertical oscillations 
around the mean orbit). These three frequencies depend on the geometry of the space-time 
and on the radius of the orbit. While they are defined as the characteristic frequencies of 
the orbital motion for a free particle, there is a direct relation between these frequencies and 
the ones of the oscillation modes of the fluid accretion flow. 

In the Kerr space-time, the radial epicyclic frequency z/^ reaches a maximum at some 
radius rmax > ''"isco and then vanishes at the ISCO, as the orbit becomes marginally radially 
unstable. Keplerian and vertical epicyclic frequencies are instead defined up to the photon 
orbit. Circular orbits with radius smaller than the one of the photon orbit do not exist. 
ug > Ur, and for corotating orbits z/r > ^9- If the space-time geometry around BH candidates 
is described by the Kerr metric, the measurement of the high-frequency QPOs may provide 
an estimate of the BH spin parameter, assuming the exact mechanism responsible for these 



oscillations is known. Torok et al. (2011) compare the available spin measurements from 



high-frequency QPOs and the continuum-fitting method and find that the two approach do 
not provide consistent results if we believe that all the observed high-frequency QPOs can 
be explained within a unique model. 

The high-frequency QPO models developed for the Kerr metric can be generalized to 
non-Kerr backgrounds (Johannsen & Psaltis 2011a Aliev et al. 2012 Bambi 2012c). The 



possibility of the existence of vertically unstable orbits, absent in the Kerr space-time, allows 
for a larger number of possible combinations between different modes, ug > Vr is not true any 
more. The conflict between measurements from high-frequency QPOs and the continuum- 
fitting method can be solved within a unique astrophysical model by admitting that the 



geometry around BH candidates deviates from the Kerr solution (Bambi 2012c). 



In four stellar-mass BH candidates, we observe two high-frequency QPOs. It is remark- 
able that in all the four cases the ratio of the two frequencies is 3:2, suggesting a strong 
correlation between them. Possible theoretical models should thus be able to explain this 
feature. The resonance models (Abramowicz & Kluzniak 2001; Abramowicz et al. 2003 



Kluzniak fc Abramowicz||2005 Torok et al.||2005" ) seem to be quite appealing from this point 



of view. For a free particle, the radial and the vertical modes are decoupled (see Eqs. (A8) 



and (A9) in Appendix). However, it is natural to expect that in a more realistic description 
there are non-linear effects coupling the two epicyclic modes. In this case, the equations can 
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be written as 



where and Fq are some functions that depend on the specific properties of the accretion 
flow, ff we knew the details of the physical mechanisms of the accretion process, we could 
write the explicit form of these two functions and solve the system. Unfortunately, that 
is not the case. However, regardless of the specific resonance model and the microphysics 
responsible for it, the resonance paradigm requires that the upper and the lower frequencies 
have the form 

vi] = miVr + m2Ve, (31) 
i/L = niVr + n2i'e ■ (32) 

where mi, m2, and n2 are integer (and likely as small as possible) numbers. Scenarios with 
a coupling between v^. and vq are surely theoretically more motivated, but it is not possible 
to exclude a priori that z/y and z/l actually arise from the combination of z/k with one of 
the other two frequencies and have the form like z^u = miz/K + m2i'r and z/l = niz/K + n2i'r- 
Assuming the JP background with deformation parameter es, we can compare the theoretical 
predictions with the observed v\] and z/l for any BH candidate and define an allowed region 
on the plane spin parameter-deformation parameter for any specific set of {mi, m2, ni, 722}- 

While systematic effects and/or wrong models are surely the most likely possibility to 
explain the disagreement between the spin measurements obtained by the continuum-fitting 



method and the high-frequency QPOs, Bambi (2012c) explores the possibility of explaining 



all the observations with a unique astrophysical model in a non-Kerr background. It turns out 
that the resonance z/g : z/j. = 3 : 1 can do the job if the JP deformation parameter is in the 
range 5 to 15. In this scenario, GRO J1655-40 and XTE J1550-564 would be "slow-rotating" 
compact objects and their ISCO would be marginally radially stable. GRS 1915-1-105 would 
instead be a "fast-rotating" object with ^ 0.5 and its ISCO would be marginally vertically 
stable. 



3.5. Jet power 



Jets and outflows are a quite common feature of accreting compact objects. In the case 
of stellar-mass BH candidates in X-ray binary systems, we observe two kinds of jets (Fender 
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et al. 2004). Steady jets occur in the hard spectral state. Transient jets appear most sig- 
nificantly when the source switches from the hard to the soft state. Most efforts so far 
concentrate on the formation mechanism of steady jets. One of the most appealing scenarios 



to explain the formation of steady jets is the Blandford-Znajek mechanism (Blandford & 



Znajek 1977), in which magnetic fields threading the BH event horizon are twisted and can 



extract the rotational energy of the spinning BH, producing an electromagnetic jet. Numer- 
ical simulations show that the mechanism can be very efficient and depends strongly on the 
BH spin ( McKinney||2005 Tchekhovskoy et al. 2010 2011). At the moment it is not clear if 
the Blandford-Znajek mechanism can be responsible for the production of steady jets, and 
in the literature there are some controversial results. The spin scenario is surely attractive, 
but still unproved. 

Recently, there have been some studies investigating if there is observational evidence for 



a correlation between BH spin and jet power in current data. Fender et al. (2010) consider 



all the spin measurements of BH binaries reported in the literature from the continuum- 
fitting method and the Ka iron line analysis. For steady jets in the hard spectral state, they 
estimate the jet power via a normalization C, defined by 

logio ^radio = C7 + 0.6 (logio " 34) , (33) 

where Lradio and Lx are, respectively, the radio and X-ray luminosity of the object. Their 
plots show no evidence for a correlation between BH spin and jet power. The left panel of 
Figure [9] shows the case in which the spin parameter is measured with the continuum-fitting 
method. If the measurements of the spins and of the jet powers are correct, it seems that 
steady jets are not powered by the BH spin. 



Narayan & McClintock (2012) propose that the Blandford-Znajek mechanism is respon- 



sible for the formation of transient jets. They consider the most recent spin measurements 
of the continuum-fitting method from the Harvard-Smithsonian CfA group, which are sup- 
posed to be more reliable, and use a different proxy for the jet power, the peak radio flux 
normalized at 5 GHz, which they claim to be model independent. These authors flnd a 
correlation between BH spin a* and jet power Pjet (see the right panel in Figure [9]), which is 
consistent with both the theoretical prediction 



Pjet OC al 



obtained in Blandford & Znajek (1977) for <^ 1 and the more accurate one 



Pjet OC OC a^/rn , 



(34) 



(35) 



where and th are, respectively, the angular frequency and the radius of the event horizon. 



found in Tchekhovskoy et al. (2010) and valid even when a^, is quite close to 1. However, nu- 



merical studies of the Blandford-Znajek mechanism show the production of steady jets, while 
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the origin of transient jets remains unclear and other scenarios may look more promising, 



such as episodic ejections of plasma blobs (Yuan et al. 2009) 



In Fender et al. (2010) and Narayan & McClintock (2012), the authors rely on spin 



measurements obtained under the assumption that BH candidates are Kerr BHs. If this 
hypothesis were wrong, all the spin measurements would be wrong too. One the contrary, 
a jet powered by the Blandford-Znajek mechanism should still be correlated with the spin 
of the compact object. The existence of the event horizon is not strictly necessary for the 
production of jets powered by the spin: for example, even neutron stars may have spin 
powered jets if there are magnetic fields anchored on the neutron star. On very general 
arguments, if the jet is powered by the rotational energy of the compact object, the power 
can be written as 



jet 



|2n 



(36) 



ra=l 



where Aq takes into account a possible non-spin contribution (Aq > 0, as Pjet cannot become 
negative for a* = 0) and the other terms are due to (some version of) the Blandford-Znajek 
mechanism. The latter must depend only on even powers of a*, because the direction of the 
spin should not be important (at least at first approximation). 

Let us now examine the spin measurements obtained with the continuum-fitting method, 
as they can be easily translated into efficiency measurements in the NT model, and then into 
allowed regions on the spin parameter-deformation parameter plane. From Section 3.2, one 



can see that an object more prolate than a Kerr BH (es > 0) may be interpreted as a Kerr BH 
with a higher value of the spin parameter. An object more oblate than a Kerr BH (63 < 0) 
has instead a spectrum similar to a Kerr BH with a lower value of a^,. While the jet power 
is independent of the spin sign, the continuum-fitting method is not, and the measurement 
of a very low radiative efficiency can be attributed to an object with a retrograde accretion 
disk. 

If we believe that steady jets are powered by the spin of the compact object, we can 
see if a correlation between jet power and spin parameter is possible for a non-vanishing 



deformation parameter (Bambi 2012d). In the left panel of Figure 9 the absence of a 



correlation between a* and Pjet is essentially due to the powerful steady jet of A0620-00, 
while the continuum-fitting method would suggest (for a Kerr metric) that this is a slow- 
rotating object. If astrophysical BH candidates were objects more prolate than Kerr BHs, 
their spin parameter would be lower than the one reported in Table [T} In particular, for 
a sufficiently large deformation parameter, A0620-00 can be consistent with a fast-rotating 
object with a retrograde accretion disk. In this case, Pjet could be high, and a correlation 
between power of steady jets and spin parameter is possible. If we neglect a possible non-spin 
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contribution to the power of the jet, a simple form for Pjet is 



Pi 



jet 



and therefore 



C = /Slogio 1*^*1 + 



(37) 



(38) 



where a' = logiQ a and /3 are the two parameters of the jet model. If we consider the JP 
metric with a universal constant common to all the objects, we find the best fit for 63 = 7.5 



a 



30.1, and /3 = 2.46 (left panel in Figure 10). If we think it is physically more plausible 



a deformation parameter that depends on the features of the compact object, the simplest 
guess is likely 



£3 = 



(39) 



with 7 constant. In this case, the best fit is for 7 = 45, a' = 31.2, and /3 = 5.65 (right panel 



in Figure 10). The puzzling result in Fender et al. (2010) may thus be interpreted as an 



indication that BH candidates are not the Kerr BHs predicted by General Relativity (Bambi 



2012d) 



If we instead believe that the Blandford-Znajek mechanism is responsible for the pro- 



duction of transient jets, as suggested in Narayan & McClintock (2012), we can constrain 



possible deviations from the Kerr geometry (Bambi 2012b). Now A0620-00 shows a quite 



weak transient jet, which is indeed what one should expect for a slow-rotating object. If we 
consider objects more and more prolate than Kerr BHs, eventually, on the base of the anal- 
ysis of the thermal spectrum of its disk, A0620-00 should be interpreted as a fast-rotating 
object with a retrograde accretion disk. This would spoil the correlation between a^, and Pjet, 
as the Blandford-Znajek mechanism would predict a powerful jet for a fast-rotating object. 

For the time being we do not know if steady or transient jets are powered by the spin 
of the compact object and therefore we cannot really use the estimate of the jet power to 
test the nature of BH candidates. If we believe that the Blandford-Znajek mechanism is 
responsible for the production of steady jets, as it is often thought, current observations 
would suggest that astrophysical BH candidates are not Kerr BHs. If we believe in the 



proposal of Narayan & McClintock (2012), in which the Blandford-Znajek mechanism is 



responsible for the formation of transient jets, despite this is not the usual scenario, we 
can put a constraint on possible deviations from the Kerr metric. The Blandford-Znajek 
mechanism can unlikely be responsible for the creation of both steady and transient jets, 
even assuming an ad hoc deformation parameter for each object, as A0620-00 shows a quite 
powerful steady jet and a weak transient jet, which would require, respectively, either a high 
or a low value of the spin parameter. 
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4. Conclusions 

Today we think tlie final product of the gravitational collapse is a Kerr BH and we 
have identified several astrophysical candidates. The latter are dark objects either in X-ray 
binary systems or in galactic nuclei, which are too heavy and compact to be, respectively, 
neutron stars or clusters of neutron stars. These two classes of objects are thought to be 
the Kerr BHs predicted by General Relativity because there is no alternative explanation in 
the framework of conventional physics. However, we have only robust measurements of the 
masses of these objects, while there is no clear indication that the geometry of the space-time 
around them is really described by the Kerr metric. 

The Kerr-nature of astrophysical BH candidates may be tested with future space-based 
gravitational wave facilities ( |Ryan|1995| [Glampedakis fc Babak|2006| [Barack fc Cutler|2007D , 



with the possible discovery of a BH binary with a pulsar companion dWex fc KopeikE|p99| 



or by studying the orbital motion of new stars at milliparsec distance from the super-massive 



BH candidate at the center of the Galaxy (Will 2008). However, all these data are not 



available now and it is definitively not clear and easy to predict when they will be available 
in the future. Recently, there have been an increasing interest in the possibility of testing the 
geometry of the space-time around astrophysical BH candidates by studying the properties 
of the electromagnetic radiation emitted by the gas in the accretion disk. These data are in 
most cases already available and, thanks to the significant progresses in the last ~ 5 years, 
people have started using these data to extract information on the space-time geometry 
around BH candidates. 

In this paper, I review the motivations, a possible approach, and some physical phe- 
nomena to probe the space-time geometry around BH candidates and check if they are the 
Kerr BHs predicted by General Relativity. The continuum-fitting method and the analysis 
of the Ka iron line are relatively mature techniques. They are now commonly used by as- 
tronomers to estimate the spin parameter of these objects under the assumption of the Kerr 
background and they can be naturally generalized to measure possible deviations from the 
Kerr solution. In general, there is a strong degeneracy between the spin and the deformation 
parameter; that is, it is not possible to measure the spin and the deformation parameter at 
the same time, but only some combination of them. That is particularly true in the case of 
the continuum-fitting method. This degeneracy can be broken (or partially solved) if we can 
constrain the space-time geometry around a specific object with two or more techniques. The 
analysis of high frequency QPOs and of the jet power are other two interesting approaches, 
but the associated physical phenomena are not yet well understood and therefore, at present, 
they cannot really be used to test the Kerr-nature of BH candidates. In the case of the high 
frequency QPOs, in a Kerr background there is no way to explain all the observations with 
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a unique scenario and find results consistent with the continuum-fitting method. While this 
fact might be interpreted as a possible indication of new physics, it is definitively possible 
that we do not have yet the correct model to describe these phenomena or that there are 
systematics effects not properly taken into account. In the case of the analysis of the jet 
power, we do not know if steady or transient jets are powered by the spin of BH candidates. 
The Blandford-Znajek mechanism can explain steady jets, but a correlation between spin 
measurement and jet power in current data would require that the space-time around BH 
candidates has deviation from the Kerr solution. If we believe that the Blandford-Znajek 
mechanism is responsible for the production of transient jets, despite this is not the usual 
scenario, observations would be consistent with the interpretation that BH candidates are 
Kerr BHs and we can put a bound on possible deviations from the Kerr geometry. 

In the near future (< 5 years), it is also expected that we will be able to observe the 
direct image of the super-massive BH candidate at the center of the Galaxy with a resolution 

. In the case of a geometrically 



thick and optically thin accretion disk, we will observe the BH shadow, a dark area over a 
bright background, whose shape depends only on the geometry of the space-time around the 



compact object and it can thus be used to test the Kerr BH hypothesis ( 


Bambi & Freese 


2009 


Bambi &; Yoshida 


2010 


Johannsen & Psaltis 


2010; 


Bambi et al. 2012 


Chen & Jing 



2012b). 
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comparable to its gravitational radius (Doeleman et al.||2008 



A. Circular orbits on the equatorial plane and ISCO 

The line element of a generic stationary and axisymmetric space-time can be written as 

ds"^ = gttdt^ + 2gtcj,dtd(j) + grr-dr"^ + gzzdz^ + g^^dcf)^ . (Al) 

Since the metric is independent of the t and (p coordinates, we have the conserved specific 
energy at infinity, E, and the conserved z-component of the specific angular momentum 
at infinity, L^. This fact allows to write the t- and 0-component of the 4- velocity of a 
test-particle as 

; _ Eg^^ + Lzgt4, : _ Egt^ + L^gu 



2 ^ r 2 

9t(f) ~ 9ttg4>4, 9t^ ~ 9ttg(i>4> 



(A2) 



From the conservation of the rest-mass, g^pX^x'^ = —1, we can write 

grr-r"^ + QzzZ^ = Ves{r, z) , (A3) 
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where the effective potential V^s is given by 

KfT = + ^^^-^^^ + - 1 . (A4) 

Circular orbits in the equatorial plane are located at the zeros and the turning points of the 
effective potential: r = i; = 0, which implies V^g = 0, and f = i; = 0, requiring respectively 
drVes = and dzVes = 0. From these conditions, one can obtain the angular velocity, E, 
and Lz of the test-particle: 



E = - ^ 9u + gt^n ^ ^^g^ 

-gtt - '^gt^^ - g.f^^'^ 



Lz = , 9t, + .,,^^ ^^7^ 

-^-gtt - '^gtcf^ - g^^^ 

where the sign + is for corotating orbits and the sign — for counterrotating ones. The orbits 
are stable under small perturbations if 9^Kff < and S^Kfi < 0. In Kerr spacetime, the 
second condition is always satisfied, so one can deduce the radius of the ISCO from S^Ks = 0. 
In general, however, one has to check also the stability along the vertical direction, as well 
as that the plunging region connects the ISCO to the surface/horizon of the compact object 



(Bambi & Barausse 2011b). 



The radial and vertical epicyclic frequencies can be quickly computed by considering 
small perturbations around circular equatorial orbits, respectively along the radial and verti- 
cal direction. If br and 5e are the small displacements around the mean orbit (i.e. r = r^^^r 
and 6* = 7r/2 + 5ei), we find they are governed by the following differential equations 



d^5e 



''+nl5r = 0, (A8) 
+ ^lSe = 0, (A9) 



df^ 
where 

The radial epicyclic frequency is Ur = Vtr/2Tx and the vertical epicyclic frequency is vg = 
VLe/2Ti. 
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log Eobs log Eobs 

Fig. 3. — Left panel: effect of the color factor /coi on the thermal spectrum of a thin disk 
for i = 45°. Right panel: effect of T on the thermal spectrum of a thin disk for i = 45° 
and 85°; IR/LDR = isotropic/limb-darkened radiation (for i = 45°, the two spectra are 
almost indistinguishable). For both panels, the parameters of the model are: M = 10 Mq, 
d = 10 kpc, M = 2- 10^^ g s-\ a, = 0.7, eg = 0.0. Flux density Ne^^^ in 7 keV"^ cm'^ s'^; 
photon energy _Eobs in keV. 



- 34 - 




0.0 0.5 1.0 0.0 0.5 1.0 

a* a. 



Fig. 4. — Xred from the comparison of the thermal spectrum of a thin accretion disk around 
a Kerr BH with spin parameter a^, and a JP BH with spin parameter a^, and deformation 
parameter e^. Left panel: = 0.7. Right panel: a,, = 0.98. See text for details. 




Fig. 5. — Constraints on the spin parameter-deformation parameter plane for the BH can- 
didates GRS 1915+105 (left panel, red area), M33 X-7 (left panel, blue area), XTE J1550- 
(left panel, yellow area), Cygnus X-1 (right panel, green area), GRO J1655-40 (right panel, 
black area), A0620-00 (right panel, magenta area). See text for details. 



- 35 - 



0.1 
0.08 
0.06 
0.04 
0.02 




a. - -6.999 

a. = 0.0 

a. = 0.7 
a. = 0.9 




1 




















a. = L 


.999 



















































2 3 4 5 6 7 
Eobs (keV) 



o 
o 

Q. 



0.1 

0.08 
0.06 
0.04 
0.02 



Eg = -5 
_ Eg = -2 
Eg = 
Eg = 2 


























£ 


g = t, 














































L 







2 3 4 5 6 7 
Eobs (keV) 



0.1 
0.08 
0.06 
0.04 
0.02 



i = 15 
- i = 30 
i = 45 
i = 60 


































1 = 75 

































































2 3 4 5 6 7 
Eobs (keV) 



o 
o 

CL 

"O 
0) 
N 

"(S 
E 
o 




2 3 4 5 6 7 
Eobs (keV) 



0.1 
0.08 
0.06 
0.04 
0.02 




i'out-i'ir = 100M 



: 1000 M 




2 3 4 5 6 7 
Eobs (keV) 



0.1 
0.08 
0.06 
0.04 
0.02 




a- = 1.313, Eg = -3 

a. = 0.837, Eg = -1 

a. = 0.7, Eg = 

a. = 0.309, Eg = 5 


















a 


• = -u 


9U 1 , i 


g - K 










































\ 







2 3 4 5 6 7 
Eobs (keV) 



Fig. 6. — Broad Ka iron line as a function of the parameters of the modeL Left top 
paneh spin parameter of the BH candidate. Right top paneh deformation parameter e^. 
Left central panel: disk's inclination angle with respect to the line of sight of the distant 
observer. Right central panel: power-law index of the intensity profile. Left bottom panel: 
size of the emission region. Right bottom panel: spin parameter and deformation parameter 
for BH candidates with the same radiative efficiency rj = 0.1036. When not specified, the 
value of the parameters is: a^: = 0.7, 63 = 0.0, i = 45°, a = —3, rout — '"in = 100 M. 
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Fig. 7. — Xrcd froni the comparison of the broad Ka iron hne generated around a Kerr BH 
with spin parameter a* and a JP BH with spin parameter a^, and deformation parameter e^. 
Left paneh a^, = 0.7. Right paneh a^, = 0.98. See text for details. 




Fig. 8. — Xredtot froKi the combination of the analysis of the broad Ka iron line and the 
continuum-fitting method for the case a^: = 0.7 (left panel) and a* = 0.98 (right panel). See 
text for details. 
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Fig. 9. — Left panel: absence of evidence for a correlation between the jet power and the 
BH spin for steady jets ( Fender et al.|[20T0 ). Right panel: evidence for a correlation between 



the jet power and the BH spin for transient jets (Narayan & McClintock 2012). See text for 
details. 
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Fig. 10. — Best fit for a possible correlation between the jet power and the spin parameter of 
BH candidates assuming a non- vanishing deformation parameter €3. Left panel: constant 
for all the objects; the best fit is for €3 = 7.5. Right panel: 63 = 71a* P, with 7 constant; the 
best fit is for 7 = 45. See Bambi (2012d) for more details. 



